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Abstract 

Using a sample of 2.59 x 10 7 ip(2S) decays collected by the CLEO-c detector, we present results 
of a study of XcO and Xc2 decays into two-meson final states. We present the world's most precise 
measurements of the X c j,(j=o,2) k + it~ , ir°ir°, K + K~ , KgK®, rjr) and 7/7/ branching fractions, 
and a search for Xc decays into 7777'. These results shed light on the mechanism of charmonium 
decays into pseudoscalar mesons. 

PACS numbers: 13.25.Gv, 14.40. Gx 
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The Xcj mesons (J = 0, 1, 2) form a triplet of cc states with one unit of orbital angular 
momentum. They are not produced directly in e + e~ annihilations, but the large branching 
fractions of ip{2S) — > Xc.il make e + e~ collisions at the ip(2S) energy a very clean environment 
for Xcj investigation. 

The Xcj mesons decay into a wide variety of different multi-hadron states. Of these, the 
two-meson states have the benefit of being comparatively straightforward to detect and to 
model theoretically. However, theoretical models based on the color singlet model make 
predictions well below the data, even when the parameters of the model are stretched to 
extremes Recent theoretical work has focussed on the Color Octet Model (COM) 
whereby contributions from the sub-process ccg — > qq are taken into account. This source 
has been shown to be a possible mechanism to make up the deficit j3|. However, these 
theoretical results were made with a very rough model for the color-octet contribution to 
the wave-function of the Xc- Fits [ij to the existing measurements [E| indicate a reasonable 
theoretical understanding of the processes involved, but data remain sparse and theoretical 
uncertainties are still large. The study of the decays to higher mass mesons (77 and 7/), offers 
the possibility of investigating the contribution of doubly-OZI suppressed decays (DOZI), 
which may compete with singly-OZI suppressed decays (SOZI) [4L and may also contribute 
to the understanding of the structure of the r\ and r/' mesons [6f] . The xa cannot decay 
into two pseudoscalar mesons because of spin-parity conservation, so we do not consider it 
further. 

The data were taken by the CLEO-c detector [7| operating at the Cornell Electron Storage 
Ring (CESR) with e + e~ collisions at a center of mass energy corresponding to the ip(2S) 
mass of 3.686 GeV/c 2 . The data corresponds to an integrated luminosity of 56.3 pb _1 and 
the total number of ip{2S) events, determined according to the method described in 0, is 
calculated as (2.59 ± 0.05) x 10 7 . 

Photons were detected using the Csl crystal calorimeter [§], which has an energy reso- 
lution of 2.2% at 1 GeV, and 5% at 100 MeV. Photon candidates were required to have a 
lateral shower shape consistent with that expected for a photon, and to be not matched with 
any charged track. We combine two photon candidates to make 7r° candidates. These are 
then kinematically constrained to the ir° mass and those combinations within 3 standard 
deviation of this mass are retained for further analysis. 

Charged particles were detected in a drift chamber system immersed in 1.0 T solenoidal 
magnetic field. The solid angle for detecting charged particles was 93% of 4tt, and the 
resolution 0.6% at 1 GeV. To discriminate charged kaons from charged pions, we combined 
specific ionizations (dE/dx) measured in the drift chamber and log-likelihoods obtained 
from the ring-imaging Cerenkov detector (RICH) [Io| to form a log-likelihood difference: 
C(K - 7r) = C RIC h(K) - C RI ch(tt) + cr 2 dE/dx (K) - cr 2 dE/dx (7c), where negative £(K - tt) 
implies the particle is more likely to be kaon than a pion. For all charged kaons we require 
C(K — 7r) < and C(K — p) < 0, and for charged pions we require C{n — K) < and 
£(7r — p) < 0. As there is potential contamination from lepton pairs in the n + ir~ final state, 
we use the muon chamber and Csl information and require that at least one of the tracks 
in this mode is inconsistent with being due to either a muon or an electron 

We use three decay modes for rj detection, 77, 7r + 7r~7r°, and 7r + 7r~7, and two modes for 
7/ detection, r/7r + 7r~ and 77r + 7r~. In each case we combine the four-momenta of the decay 
products into an 77 W candidate, kinematically constrain the candidate to its nominal mass 
and retain those candidates with a fit x 2 < 9/1 degree of freedom. 

For events with two distinct meson candidates, we combine the candidates into a Xc 
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candidate. At this stage of the analysis, the invariant mass resolution of the Xc is approx- 
imately 15 MeV/c 2 . We then search for any unused photon in the event and add that to 
the Xc candidate to form a ip{2S) candidate. This if)(2S) candidate is then kinematically 
constrained to the four-momentum of the beam, the energy of which is calculated using the 
known ip{2S) mass. The momentum is non-zero due to the finite crossing angle (~ 3 mrad 
per beam) in CESR. To make our final selection, we require the ip{2S) candidate to have a 
X 2 of less than 25 for the four degrees of freedom for this fit; this requirement rejects most 
background combinations. This kinematic fit greatly improves the mass resolution of the Xc 
candidate to values ranging from 3.2 to 8.7 MeV/c 2 depending on the spin of the Xc and the 
decay mode. 

To study the efficiency and resolutions, we generated Monte Carlo samples for each Xc into 



each final state using a GEANT-based detector simulation The simulated events have 
an angular distribution of (1 + a cos 2 8), where 6 is the radiated photon angle relative to the 
positron beam direction, and a = 1 and 1/13 for the Xco, and Xc2 respectively, in accordance 
with expectations for an El transition. The decay products of the Xco are generated with 
a flat angular distribution. The products of the Xc2 are generated according to a double 
correlation function of the polar angles of the mesons measured in the Xc rest frame relative 
to the transition photon direction [12j. The efficiencies are shown in Table I. The efficiencies 
for i] and r]' modes include the relevant branching fractions. 



TABLE I: Yields found in the data sample and detection efficiencies obtained from analyses of 
Monte Carlo generated events. 



Mode 


XcO 


Xc2 




Yield 


Efficiency(%) 


Yield 


Efficiency(%) 




8934 ± 111 


58.7 ±2.4 


2543 ± 56 


66.2 ±2.7 




2807 ± 62 


40.0 ±4.4 


793 ± 33 


48.5 ±5.3 


K + K~ 


8156 ± 100 


53.8 ±2.5 


1645 ± 42 


60.2 ±2.8 




2109 ± 49 


25.3 ± 1.1 


373 ± 20 


29.3 ± 1.3 


r]rj 


930 ± 35 


12.3 ±1.1 


156 ± 14 


12.6 ± 1.1 


r/rj' 


35 ± 13 


9.2 ±0.8 


3.3 ±8.0 


10.5 ±0.9 


rfrf 


413 ± 24 


8.2 ±0.6 


12 ±7 


8.8 ±0.5 



The final invariant mass distributions are shown in Fig. 1. These distributions are then 
fit with two signal shapes comprising Breit-Wigner functions convolved with Gaussian res- 
olutions, together with a constant background term. The masses and widths of the Breit- 
Wigner functions were fixed according to the Particle Data Group (PDG) averages and 
the widths of the Gaussian resolution functions were fixed at the values found from Monte 
Carlo simulation. The yields from these binned likelihood fits are tabulated in Table I. 

To convert the yields to branching fractions, we divide by the product of the number 
of ip(2S) events in the data sample, the detector efficiency, and the branching fractions for 
ip(2S) into Xcj- For the last factor we use the CLEO measurements of BU])(2S) — > 7Xco) = 
(9.22 ± 0.11 ± 0.46)% and B(i/>(2S) -> 7^2) = (9.33 ± 0.14 ± 0.61)% Q. The results are 
tabulated in Table II. 

We consider systematic uncertainties from many different sources. All modes have a 
2% uncertainty from the total number of ip{2S) decays |§]. The requirement on the x 2 °f 
the constraint to the beam four-momentum has been checked by changing the cut value 
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in the range 12-50 and noting the change in the yield in these, and other similar decay 
modes. Based on this study we place a systematic uncertainty of 2.5% on the efficiency of 
this requirement. The uncertainties due to track reconstruction are 0.3% per charged track 
(0.67% for kaons). The limited Monte Carlo statistics introduces an uncertainty which is 
in all cases less than 1.5%. The systematic uncertainty due to the photon detection and 
shower-shape criteria is set at 2% per photon both for the transition photon and for the 
decay products of rj and tt° decays. In the cases including 77 decays, this contribution is 
incorporated taking into account the fraction of those decays that proceed through each rj 
decay mode. The final signal plots are all well fit using the functions described above. By 
studying the variation of the yields of the high statistics modes resulting from floating the 
signal parameters, we assign a 2% uncertainty in each mode due to the uncertainties in the 
fitting procedure. We have checked that the yields from the various decay modes of the 
77 W mesons are consistent with their branching fractions and efficiencies. When calculating 
the final branching fractions, we add the above systematic uncertainties in quadrature. The 
uncertainty due to the ip(2S) — > 7Xc branching fractions is kept separate and quoted as a 
second systematic uncertainty. 

For evaluating the limits in the cases where there is no significant signal, we take the 
probability density function and convolve this with Gaussian systematic uncertainties. We 
then find the branching fraction that includes 90% of the total area. 



TABLE II: Branching fraction results (in units of 10~ 3 ) for each decay mode. The uncertainties 
are statistical, systematic due to this measurement, and systematic due to the ip(2S) — ► XcJ7 rate, 
respectively. The limits on the branching fractions include all systematic uncertainties, and central 
values for those measurements are included in parentheses. 



Mode Xco Xc2 



7T + 7T 


This Work 


6.37 ±0.08 ±0.29 ±0.32 


1.59 ±0.04 ±0.07 ±0.10 




PDG [5] 


4.87 ± 0.40 


1.42 ±0.16 


TT°TT° 


This Work 


2.94 ±0.07 ±0.32 ±0.15 


0.68 ±0.03 ±0.07 ±0.04 




PDG 


2.43 ± .20 


0.71 ±0.08 


K+K- 


This Work 


6.47 ±0.08 ±0.33 ±0.32 


1.13 ±0.03 ±0.06 ±0.07 




PDG 


5.5 ±0.6 


0.78 ±0.14 




This Work 


3.49 ±0.08 ±0.17 ±0.17 


0.53 ±0.03 ±0.03 ±0.03 




PDG 


2.77 ±0.34 


0.68 ±0.11 


T]7] 


This Work 


3.18 ±0.13 ±0.31 ±0.16 


0.51 ±0.05 ±0.05 ±0.03 




PDG 


2.4 ±0.4 


< 0.5 


f]f]' 


This Work 


< 0.25 

(0.16 ±0.06 ±0.01 ±0.01) 


< 0.06 

(0.013 ± 0.031 ± 0.001 ±0.0( 




PDG 


< 0.5 


< 0.26 


rj'rj' 


This Work 


2.12 ±0.13 ±0.18 ±0.11 


< 0.10 



(0.056 ± 0.032 ± 0.005 ± 0.003) 
PDG 1.7 ±0.4 <0.4 



fromBES 14 and CLEO 15 



to results 
2B{ X c - 

Our results do not include that constraint, but the data are consistent with this 
isospin symmetry. Our results are also consistent with the expected result that B(xc 



A 



Our results are summarized in Table II, and compared with the PDG fits |5J 

These fits explicitly assume that B(xc — ► 7r + 7r~) 
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K + K~) = 2B(xc -^5-^5)5 whereas previous results had indicated that this might not be 
so in the J = 2 case. The largest deviation from previous results (~ 3er) is in the case of 
Xco ~^ 7i + 7i~. In the case of the Xc2, our limit for the branching fraction into 777/ is below the 
fit value obtained from previous data by Qiang Zhao [If, suggesting that the DOZI decays 
of the Xc2 m ay contribute less than indicated by that phenomenological analysis. We note 
that there is an overlap of datasets in the results presented here and those of our previous 
analysis of rf'^rf- 1 ' decays, and so our new results should replace rather than augment our 
previous measurements. 

In summary, we measure branching fractions for Xco and Xc2 decays into 7r°7r°, 7r + 7r~, 
K + K~, K^Kg, r]T], and rfvj '. The decay Xc2 W1 * s observed for the first time and in all 
other cases these measurements are more precise than any previously made. These results 
may be used to test the role of the Color Octet Mechanism model of Xcj decays. The 
improved limits on decays into r/r/' are further proof of the small contributions made by 
DOZI decays in this system. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, the Natural Sciences and Engineering Research Council of Canada, and 
the U.K. Science and Technology Facilities Council. 
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FIG. 1: Invariant mass distributions for ir + ir , 7r°7r°, K + K 



K° S K° S , r/r/, 7/7/, 7/7/ 



The fits are 



described in the text. The downward arrows are at the value of the invariant mass of the Xci- 
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